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(57) Abstract 

We present a method and apparatus for local 
and superficial measurement of the optical proper- 
ties of turbid media. The depth probed is on the or- 
der of 1 transport mean free path of the photon. The 
absorption coefficient, rp.fiiirp.rl scatte ring c.qeS^c\ejvX 
And^the^phase^function^parameter^ 7=( i-g2)7(T~-g7F 
are optical parameters computed from a single mea- 
surement of the spatially resolved reflectance close 
to the source. Images of superficial structures of the 
medium can be obtained by performing multi-site 
measurements. An important application of this 
technique is the characterization of biological tis- 
sues, for example for medical diagnostic purposes. 
Measurements on biological tissues can be achieved 
using a probe of diameter less than 2mm, and the 
average volume probed is on the order of 1 mm 3 . 
Separate images of absorption and tissue structure 
can be achieved with a resolution of approximately 
one transport mean free path of the considered tis- 
sue. 
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METHOD AND APPARATUS FOR MEASURING LOCALLY AND SUPERFI- 
CIALLY THE SCATTERING AND ABSORPTION PROPERTIES OF TURBID 
MEDIA 



The present invention relates to a method and an apparatus to quantify the optical scatter- 
ing and absorption properties of a turbid media. More precisely the present invention 
relates to a non-invasive measurement, over a small area of the sample surface. Local and 
superficial characterization of biological tissues in vivo is a major application of this 
invention. 

Z Related Background Art. 

Different techniques have already been proposed to quantitatively determine the absorp- 
tion and reduced scattering coefficients of turbid media 1 . Most of the non-invasive meth- 
ods are based on the measurement of spatially and/or temporally-resolved reflectance. The 
principle is as follows: the turbid medium is illuminated by a light source. The backscat- 
tered light is measured by one or several detectors. Different types of measurements are 
possible, depending on the time-dependence of the illuminating source: steady-state (con- 
tinuous source), time-domain (short pulsed source) or frequency domain (amplitude mod- 
ulated source). The present invention relates to the case of steady-state measurements, 
performed at different distance p between the source and the detectors. However, the tech- 
nique presented here can be complemented by time- or frequency-domain measurements. 
The range of p values is an important point to consider, when comparing different meth- 
ods based on the measurement of the reflectance. First, the probed volume of the turbid 
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medium is related to the source-detector separation p. The larger the source-detector sepa- 
ration, the deeper the average depth probed. Second, depending on the range of p, different 
mathematical processing must be used to obtain the optical properties from the raw data. 
At least two cases must be distinguished. 
^ 1) The first case corresponds to source-detector separations larger than several transport 
mean free paths. For typical biological tissue optical properties 2 , this case correspond to 
source-detector separation of larger than 2 mm. An analytical form of the reflectance can 
be obtained from the diffusion equation, if the absorption coefficient p a is sufficiently 
lower than the scattering coefficient |a s (typically at least ten times). In such a case, the rel- 
10 evant optical properties are the refractive index, the absorption coefficient and the reduced 
scattering coefficient. The average depth of probing is on the same orde p than the 
source-detector separation. 

Such methods have been already published, and are the object of patents (Ref.l, Patent 
5,517,987 Tsuchiya, Patent 5,676,142 Miwa et aL). 
15 2) The second case corresponds to source-detector separations close to one transport mean 
free path. For biological tissues 2 , such source-detector separations correspond typically to 
distances from 0.1 to 2 mm. The average depth probed is on the order of 1 mm. Such small 
source-detector separations enable the measurement of the optical properties of a small 
tissue volume. 

20 In this case, light propagation can be modeled using Monte Carlo simulations. In contrast 
to the previous case 1), it was discovered by the inventors that not only the first moment of 
the phase function must be considered, but also the second moment. More precisely, the 
relevant optical properties are the refractive index, the absorption, the reduced scattering 
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coefficient and the parameter y=(l -go V(l-giX that takes into account the first two moments 
of the phase function, denoted g { andg 2 respectively. 

Different methods have been proposed for local characterization of turbid media, and in 
particular of biological tissues. Nevertheless none allows for the simultaneous determina- 
5 tion of the absorption, the reduced scattering coefficient and the parameter y. Wang et al. 
(Patent 5,630,423) proposed a method for the determination of the reduced scattering 
coefficient only, using an optical beam of oblique incidence. Moreover, their analysis does 
not include the effect of the phase function. Kessler et al. (Patents 5,284,137 and 
5,645,0619) proposed a method for measuring the local dye concentration and scattering 
10 parameters in animal and human tissues, based on spatial and spectral measurements. 
However, their methods do not enable the simultaneous determination of the absorption 
coefficient, reduced scattering coefficient and the parameter y. 

In the present invention we present a method and apparatus for the measurement of the 
absorption coefficient, the reduced scattering coefficient and the said phase function 
15 parameter y, from the spatially-resolved reflectance data at short source detector separa- 
tion. These three parameters, which can be measured at different wavelengths, enable us to 
characterize turbid media, such as biological tissues. Our method is based on a theoretical 
study of the reflectance at short source-detector separations, that we performed with 
Monte Carlo simulations. Other publications may be of concern: 

20 1. See "Welch, A. J.; van Gemert, M. J. C. Optical Thermal Response of Laser Irradiated 
Tissue; Plenum publishing Corp., New York, 1995", and references therein. 

2. W.-F. Cheong, S.A. Prahl, and A.J. Welch, "A Review of the Optical Properties of 
Biological Tissues," IEEE J. Quantum Electron. 26, 2166-2185 (1990). 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1. Description of the principle of spatially-resolved reflectance measurement. 

Fig. 2. Probability density function of the mean depth of scattering events per detected 
photons. 

Fig.3. Examples of reflectance obtained with different phase functions. Case of matched 
D refractive index (n=1.0). 

Fig.4. Monte Carlo simulations and fits on the form /?(p) = [A(p, \x s ; y) + ix s 'B(ix 0 )f . Case of 

constant y=l .9 and fi s '=l mm" 1 , numerical aperture of the source and 

detectors = 0.37. A = 0.0647p°* 324 exp(-0.16lp), B = 0.18653 |a a - 0.8466 |u a 2 + 1.836|i a 3 

Fig. 5. a. Basic description of the apparatus: Case of measurement with an optical fiber 
probe 

10 

Fig. 5.b. Basic description of the apparatus: Particular embodyment where the probe is 
composed of n detecting fibers (typically 6). A broadband or white light source is used for 
excitation and the retrodiffused light is collected by the n detecting fibers and dispersed in 
parallel in a spectrograph, before being detected on a 2D CCD camera and transmitted 
numerically to a PC for y, \i s \ jn a determination at each wavength X . 

j ^ Fig. 5.c. Basic description of the apparatus. Case of measurement with sources and 
detectors directly in contact with the turbid medium. 

Fig. 5. d. Basic description of the apparatus. Case of non-contact measurements with a ID 
or 2D detector, coupled with an imaging device. Scanning systems 1 and 2 are optional, 
they may be distinct from one another or confounded in a single scanner. 

Fig. 5.e Basic description of the apparatus. Case of non-contact measurements with 
20 source and detection beams deflected by one or two scanning devices which may be 
distinct from one another or confounded. 

Fig. 6. Examples of the sample side of the fiber optical probe. 

6a. simple arrangement for a single measurement. 
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6b. arrangement for symmetrical measurements, 
6c. Arrangement for multiple measurements. 

6d. Arrangement for multiple measurement, using an multicore fiber, or an optical fiber 
bundle. 

5 6e. Arrangement for symmetrical measurements where detecting fibers are arranged on an 
elliptic path so as giving a regular spacing from the illuminating fibers to the detecting 
fibers, while maximizing the collecting surface of the fibers. Additional fibers F can be put 
in the center of the probe or around the probe for other light channels: collecting fibers for 
fluorescence or Raman scattering for example. 

^ Fig.7.ExampIe of an measured the spatially-resolved reflectance with the probe shown in 
Fig. 6a. The measurement, performed on a microsphere suspension, is superimposed to a 
Monte Carlo simulation. The optical properties, computed from published water optical 
properties (ji a ) and Mie theory (|n s \ y) are: n=1.33, ^=0.0004 mm" 1 , |V=1.0 mm' 1 , y=2.2. 
The calibration was performed with a siloxane sample of known optical properties. 

Fig.8.Relation between the parameters R(p=l mm) and |9plnR(p=l mm)| and the optical 
15 coefficients |a s ' and n a Case of y=1.5 and 1.9. Probe of refractive index = 1.5, sample of 
refractive index = 1.4, optical fibers diameter = 200 \im, NA=0.37 (source and collection). 

Fig.9.Plot of the parameter JLjR for different y values (1.0, 1.8, 2.5), different reduced 
albedo a' (1, 0.95, 0.9) and for fixedp=l mm. Mismatched refractive index n=1.4. 

Fig.lOPlot of jR(a= \)-jR( a y for different y values and reduced albedo a\ and for fixed 
p=l mm. Mismatched refractive index n=1.4. 

20 
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DETAILED DESCRIPTION OF THE INVENTION 

I. Definitions 

The concept of spatially resolved reflectance is illustrated in Fig.l. Consider light imping- 
ing on a turbid medium (through air or through a light guide), in a given solid angle 
^source- The spatially resolved reflectance R(p) is defined by the backscattered light power 
in a given solid angle cOdetecton at a distance p from the source, per unit area and normal- 
ized by the source power. The distance p is referred as the source-detector separation (also 
in the case when light guides or imaging devices are used for the illumination and for the 
collection of the backscattered light). 

R(p) depends on the optical properties of the turbid medium, defined below. Note that 
R(p) depends also on the source and detectors characteristics, i.e numerical aperture and 
sizes. 

It is commonly admitted that the fundamental optical properties of a turbid medium is 
determined by the average index of refraction n of the medium, the absorption coefficient 
|i a , the scattering coefficient n s , and the phase function p(9) where 9 is the scattering 
angle. The absorption coefficient ja a [mm" 1 ] is defined as the probability of absorption per 
unit infinitesimal pathlength. The scattering coefficient (a s [mm" 1 ] is the scattering proba- 
bility per unit infinitesimal pathlength. The phase function p(9) is the density probability 
function for the scattering angle 9. The phase function is normalized as follows: 



20 1 = 2n\ p(Q)s\nQclB . (1) 



The n th order moment g n of the phase function is defined as: 

gn = 27if n p„(e)p(e)sin0^e (2) 
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where P n (0) is the Legendre polynomial ot order n. The first moment of the phase function 
is also called the anisotropy factor, and is often simply noted g (=gj). It represents the 
mean cosine of the scattering angle. The reduced scattering coefficient (a s ? is defined as: 

M s ' = Hs(l-gj) (3) 
5 The transport mean free path (or reduced mean free path) mfp' is defined as: 

mfp^OV+Ma)" 1 (4) 

The reduced albedo a' is the ratio: 
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a , =ii s 70i s '+ n a ) (5) 
As a result of the present invention described in more details in the next section, it is also 
necessary to define another phase function parameter called y: 

Y=(l-g 2 )/d-gi) (6) 
All the parameters listed above are referred as optical properties. They are wavelength 
dependent, and can vary in space and time. 

2i Reflectance measurements at distances close to one transport mean free path 
The methods described in the Patent 5,517,987 (Tsuchiya) are based on measurements 

with large range of source-detector separations, typically from 1 mfp' to 10 mfp\ In such 

cases, the volume probed is on the order of 10-1000 (mfp') 3 . In contrast with such a large 

scale investigation, the present invention is directed to a novel approach where the volume 

probed is much smaller, on the order of 1 (mfp') 3 . This is achieved by using only small 

source detector separations, typically from 0.1 to 2 mfp\ The lateral dimension of probing 

is limited to this range of distances. 
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Note that the present invention can be complemented by large source-detector separation 
measurements, if the optical properties of both superficial and deep pans of a turbid 
medium must be determined. 

A model of photon migration in tissue was necessary to predict the relationship between 
2 the measured reflectance and the optical properties. Analytical solutions from the diffusion 
equation are not appropriate in our case because we are interested in the reflectance close 
to the source, at distances comparable to the transport mean free path [mfp']. We have per- 
formed Monte Carlo simulations to predict the measured reflectance of an homogeneous 
semi-infinite turbid media. 

10 The exact diameter of the illuminating and collecting fibers, as well as their numerical 
apertures, have been taken into account in the simulations. The mismatch of index of 
refraction at the surface of the medium have been also taken into account in the simula- 
tions, by using the Fresnel law for each photon reaching the surface. 
This is also one result of this simulation to compute the average depth of probing, illus- 

15 trated in Fig.2. It is demonstrated that only the superficial part of the turbid medium is 
probed if small source-detector separations are used. For this, we determined the depth 
below the surface of each scattering event in the simulation. With this information we 
determined the average depth of all the scattering events for each detected photon, which 
we present as a probability density function in Fig.2. This figure shows that the average 

20 depth of scattering is approximately around 1 mfp'. Moreover it showed that the part 
located below 2 mfp' were not playing a significant role in the measured signal (for p< 1.5 
mfp'). 

The spatially resolved reflectance R(p), with short source-detector separations is more 
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complex than in the case of large source-detector separations. Indeed, for small source 
detector-separation conditions, the inverse problem, i.e. calculating the localized absorp- 
tion and reduced scattering coefficients, is sensitive to the scattering phase function. It is 
pan of the present invention to have shown, from Monte Carlo simulations, that only the 
2 first and second moments of the phase function can be taken into account. Moreover, it 
was established that the influence of the these two moments are not independent. Indeed, it 
is a merit of this invention to demonstrate that only one parameter Y=(l-g2y(l-gi)» which 
depends on the first and second moments of the phase function can be used to characterize 
accurately the reflectance at short source-detector separation. 

10 Fig.3 illustrates the fact that the parameter y is the only predominant parameter of the 
phase function that must be taken into account (and not g t as frequently mentioned in lit- 
erature). Different reflectance curves, obtained from Monte Carlo simulations are shown 
in Fig.3. Four different phase functions were used for the simulations. Three phase func- 
tions are characterized by y=1.25, but different gj and g 2 values. Fig.3 shows that almost 

15 identical reflectance curve are obtained for distances p|i s '>0.3, if the gj and g 2 are varied 
in such a way that the parameter y stays constant (y=1.25). For comparison, the reflectance 
computed with a different g value is also presented (y=2.25). It appears then clearly that 
the parameter g plays a significant role in the reflectance. 

The parameter y depends on the characteristics of the set of scatterers inside the turbid 
20 medium: shape, distribution of sizes and refractive index. 

The possibility of determining this parameter y is an important part of this invention, and 
has never been reported earlier. The determination of y can be related to the microscopic 
structure of the sample, since y is related to the sizes and refractive indexes of the scatters. 
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In most cases, the average refractive indices n and the source and detectors characteristics 
are fixed and known. In such case, only three parameters are determined from measure- 
ments at short source-detector separation: ^t a , |i s ' and y. 

It is important to note that the parameter y, cannot be estimated if large source detector 
separations are used. Indeed, the reflectance is only sensitive to y at short source-detector 
separation. 

It was also derived from Monte Carlo simulations that the reflectance R(p), for p=lmfp\ 
can be reasonably well approximated by the expression: 

R(9) = M(p, \x s \ y) + B{\i a , [i/)f (7) 

In many cases, we found that Equ.(7) can also be written as: 

*<p) = [A(p, n;, 7 ) + v s 'B(ix a )f (8) 
it is important for the description of the invention to evidence the fact that the function A 
depends on p, the scattering properties (i.e. \i s * and y) but not on |Li a . In contrast, the func- 
tion B depends on ^ a and ji s ' but neither on y nor p. An example of Equ.(8) is shown in 
Fig.4. 

This formulation is of great help to solve the inverse problem, as described in section 5.3. 
3,. Apparatus 

In the first embodiment, The apparatus can be divided in three parts, described in Fig. 5. a. 
The first part is the illuminating system. Any light source can be used. For examples: a) 
white sources, such as halogen or xenon lights, metal halides or fluorescent or phosphores- 
cent sources, b) sources such as lasers, laser diodes, optical fiber lasers, light emitting 
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diodes or superluminescent diodes, c) the sources described in point a) and b) where 
monochromators, filters or interference filters are added to select a given set of wave- 
lengths. 

In the first preferred embodiment, The light power is conducted to the investigated sample 
5 by the probe, which is the second part of the apparatus. The probe is preferably made of 
optical fibers, to illuminate and to collect the backscattered light. But GRIN rods or other 
type of light pipes can also be used. Different possible arrangements of optical fibers are 
illustrated in Fig.6. Two different modes of measurements can be chosen. First, one fiber is 
used to illuminate the sample and at least two other fibers are used to collect the backscat- 
10 tered light at two different distances. Second, one fiber is used to collect the light and at 
least two other fibers, located at two different distances from the first one, are used to illu- 
minate sequentially the sample. 

The arrangement of the different fibers can be replaced by any imaging system or image 
guide, such as multicore optical fibers, or optical fiber bundle. 

15 The light collected by the probe is analyzed by the detection unit, which is the third part of 
the apparatus. If wide spectral light sources are used (such as halogen or xenon lights), a 
spectrograph can be put between the probe and the detector to get wavelength dependence 
of the backscattered signal (either in the source or detection unit). Any types of detectors 
can be used. For example, photodiodes, avalanche photodiodes or photomultipliers can be 

20 assigned to each collecting fibers. Simultaneous detection of each collecting fiber can also 
be achieved using linear or two-dimensional detectors such as Charge-Coupled Detectors 
(ID or 2D), intensified CCD or array of photodiodes. 

Aparticular embodiment is described in Fig. 5.b. : the probe is composed of n detecting 
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fibers (typically 6). A broadband or white light source is used for excitation and the retrod- 
iffused light is collected by the n detecting fibers and dispersed in parallel in a spec- 
trograph, before being detected by a 2D CCD camera and transmitted numerically to a PC 
for y, H s \ Ma determination at each wavelength X.This particular embodiment has the 
5 advantage of yielding the wavelength dependance of y, |i s \ M a frorn a single measurement. 
A second type of embodiment is presented on Fig. 5x. The difference with embodiments 
presented in Fig. 5. a. and Fig. 5.b. is that optionally no optical fibers, light pipes or grin 
rods is used. The light source unit is directly in contact with the turbid medium, as well as 
the detector unit. Collimating optics, microoptics or imaging optics (DOE for example) 

10 can be put between the turbid medium and the actual light sources and detectors. The dif- 
ferent type of sources and detectors cited in example for the first embodiment can be used 
for the second type of embodiments. Hybrid design, such as arrangements involving both 
direct contacts sensors or detectors and fibers, light pipes or grin rods are also included in 
the present embodiment. 

15 The third embodiment is described in Fig. S.d.and Fig. 5.e. The only difference with the 
first and second embodiment is that no contact measurements are performed. A collimat- 
ing. system allow for a point-like illumination on the turbid medium. An imaging system 
enables the measurement of the spatial distribution of the reflectance. The detectors can be 
either an array (ID or 2D) of detectors (Fig. 5.d) with optional scanning 1 and 2 which can 

20 operate separately or can be confounded in a single scanner, or a single detector (Fig. 5.e). 
In this last case, an scanning device is used to obtain the spatially-resolved reflectance. A 
fiber bundle, multicore fiber or relay optics (grin rod or multiple lenses) can be put 
between the focal point of the imaging system and the detector(s). The different type of 
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sources and detectors^cited in example for the first embodiment can be used for the third 
embodiment. 

Spatial images of the parameters |i a , |u s ' and y can be obtained by a series of measurements 
at different locations, that we call multi-site measurements. Some mechanical or optical 
- scanning device can be used for this purpose. The resolution of such images is on the order 
of the mean source-detector separation used for a single site measurement. All three 
embodiments can be expanded to perform multi-site measurements, by duplicating and/or 
multiplexing the illuminating or measuring devices. For example, the optical fiber probes 
shown in Fig. 6, can be duplicated and put side by side. The scanning system described in 
jq the third embodiment can also be expanded to perform multi-site measurements. 



4. Normalization and calibration 

The differences of transmission between each fiber are corrected by performing a mea- 
surement on a turbid phantom illuminated uniformly. 

The background light, measured with the light source turned off, must be subtracted from 
the signal. 

In order to perform absolute intensity measurements, a calibration is performed on a turbid 
medium of known optical properties. Examples of such a medium are: 1) solid or liquid 
turbid medium which properties have been measured by other techniques, or reported in 
the literature 2) water suspension of microsphere of known size distribution and refractive 
index. Absorbing dye may be added to the suspension. In case 2) the scattering properties 
are calculated from Mie theory, and the absorption coefficient is assumed to be equal to 
the water absorption coefficient, if no absorbing dye is added. If an absorbing dye is used, 
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the absorption coefficient can be measured by a spectrophotometer, before mixing the 
solution with any scattering materials. 

A Monte Carlo simulation is performed with the optical properties of the calibration sam- 
ple. The simulation is then divided by the experimental reflectance performed on the cah- 
^ bration sample. The result, that must be independent of the source-detector separation, is 
defined as the calibration factor. Each new measurement is multiplied by the calibration 
factor. 

An example of a measurement, after calibration, is shown in Fig.7. The measurement was 
performed with the apparatus described in the first embodiment, with a probe similar to 

10 the one described in Fig. 6. a. A prior calibration was performed on a solid turbid medium, 
which optical properties were measured by another technique (frequency domain mea- 
surement, cf. Refl). The sample was a water suspension of polystyrene microsphere of 1 
jam diameter. The measurement is compared to a Monte Carlo simulation performed with 
the scattering properties calculated from Mie theory and absorption coefficient of water. 

15 Fig.7 shows an excellent agreement between the measurement and the simulation, which 
confirm the validity of our simulation model, experimental measurement and calibration. 

5. Signal processing 
5.1. Control of the homogeneity of the area probed 
20 Artifacts during a measurement, for example due to bad contact between the probe and the 

sample, or heterogeneity of the sample, can be detected by the following optional proce- 
dure. Two illuminating fibers are disposed symmetrically in regard to the collecting fibers 
(see Fig.6b). If the sample is homogeneous, the reflectance curve should be identical with 
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either illuminating fiber. Therefore, heterogeneity of the investigated region or obstruc- 
tions beneath the fibers are detected by comparing the two curves. If the two curves are 
sufficiently close, the measurement is validated and the average of the two curves is calcu- 
lated. 

5.2. Smoothing procedure and computation of the derivative of the curve . 
Functions in the form m { p m2 exp(m 3 p) were always found to fit well Monte carlo simulations 

for a restricted range of distances. Therefore, a smoothing of the experimental reflectance 

R(p) can be processed by fitting R(p) to mjp^exp^p) . 

The determination of the slope of the logarithm ~-{\nR{Q, \x/ t a\ y)) is also derived from this 
fit, using the following formula: 



The slope of the square root of R(p) can also be derived (the use of these slopes are 
described below in section 5.3): 



5.3. Inverse problem 

We developed different methods to solve the inverse problem, which consists in extracting 
optical coefficients from the reflectance data. 



Monte Carlo simulations show that the measurements of the reflectance intensity R(p) and 
the slope of lnR(p) (denoted 3plnR), determined at a fixed distance p, can be used to 
derive |a s ' and for a given y value. Fig. 8 shows graphically the relationship between ji s ' 




(9) 




(10) 



Method 1. 
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and |i a and the two parameters R(p=l mm) and |9plnR(p=l mm)| in the case of for 
p=l mm, y = 1.5 and 1.9. Note that any other choice of p is possible, as long as p is on the 
order of I mfp'. We see in Fig.8 that ja s ? and ji a can not be determined uniquely if y is 
unknown. Further insight for optimizing the inversion strategy is provided by three addi- 
5 tional features of Fig.8. 

First, the determination of \i s y is only weakly influenced by y. Indeed in Fig.8 the errors 
induced by error in y are typically 10% for |i s ' around 1 mm' 1 . Second, although absolute 
determination of is not possible when y is not precisely known, relative absorption 
changes can be still evaluated. Third, the indetermination of the parameter y may be 
10 resolved by the values of R(p) and/or |8plnR(p)| at other distances. Therefore the follow- 
ing procedure can be used: 

(1) determination of |a s ' and ja a from R(p=l mm) and |dplnR(p=l mm)| for a set of values 
Y- For example: y= 1.0, 1.1, 1.2,. ..,2.5. 

(2) simulations with the different sets of |i s ' and ji a obtained from point 1) 

15 (3)comparison between the simulations and the reflectance profile for distances 0.3<p<2 
mm. 

This last step allows us to determine the correct values of y, \i s * and ju a . Points 1 to 3 can 
be done iteratively to evaluate y more precisely, using a finer discrimination of y values. 
This method described by the iteration of points 1) to 3) corresponds to finding a simula- 
20 tion that fit closely a measured curve R(p). This can also be achieved by directly compar- 
ing the measured curve R(p), with a set of simulations, and finding the simulation that 
minimize the square of the differences between the measured and simulated curves. 
Note also that the use of the reflectance intensity R(p) and the slope of lnR(p) is equivalent 
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to the use of two intensities R(Pi) and R(p2>, measured at two close distances p x and p^, 
respectively. 

Method 2. 

5 The inverse problem can also be performed, considering properties of R(p) expressed by 
Equ.(7). 

Indeed, it can be derived from Equ.(7) that the quantity JLy*(p, n/, \x a .y) does not depend 
on the absorption coefficient |i a : 

^7/?(p,M/,^Y)= 2§£(p.n 5 ', Y ) (11) 
^ This property is confirmed in Fig.9, where the quantity ^JR(P, h/,y), evaluated at p=l 
mm, is plotted as a function of fi s * for y =L 1.9 and 2.5 and reduced albedo a'=l, 0.95 and 
0.9 (note that any other choice of p is possible, as long as p is of the order of 1 mfp f ). 
In Fig.9, the parameter Ay^p, \x/ t y) clearly depends on \i s * and y. In contrast, the depen- 
dence on a' is almost negligible. 
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Therefore, y and (i s ' can be derived from the parameter JLjR{p) , calculated from the 

dp 

experimental reflectance R(p) (see section 5.2). Simultaneous determination of y and |i s ' 
require values of j^-Mp) at different distancesp (at least two). If one of the parameters y or 
|i s ' is already known, the determination of the other one can be obtained from the value of 



3pV*(p) at a single distance. 



Convenient analytical approximation of ^V*(P) can be obtained by fitting Monte Carlo 



results to polynomial functions. For example, in the case of fixed y=1-9, we obtained: 

H s '=0.9162-52.89x+1795x 2 -18155x 3 +65428x 4 (12) 
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where x= \jJL<Mi>)) at p=l mm. x is expressed in [mm" 2 ] and (J s y is expressed in [mm* 1 ]. 

This results is valid for a probe of refractive index of 1.5 and a sample of refractive index 

of 1.4, optical fibers of diameter 200 pm, NA=0.37 (source and collection). 

Once p s ' and y are calculated from the procedure explained above, |j a is calculated from 

the absolute value of R(p), which highly depends on For given \i s \ y and p values, the 

dependence of reflectance on ja a is obtained by Monte Carlo results. From Equ.(7) we 

have: 

n a = mVrcpI-Ay,^)] (13) 
where / and h are functions given by Monte Carlo simulation. Particularly, they can be 
well approximated by polynomial functions. For example, for y=\.9, probe refractive 
index of 1.5, sample refractive index of 1.4, optical fibers of diameter 200 jam, NA=0.37 
(source and collection): 

h =-0.002257-8. 171y+268.8y 2 (14) 
/= 0.01311+0.05184^ s , -0.01974 t i s ' 2 + (15) 
0.003217^l s ' 3 -0.0001992|Li s ,4 



where v = [-/R(p)-/] [mm -1 ]. ja a is expressed in [mm" 1 ]. 



Method 3 

20 

Equ.(7) also shows that relative measurements of n a , i.e. variation of ja a from a known 
value ^i ao , is possible by measuring the variation of the parameter jR{p) : 
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This relation is illustrated in Fig. 10 in the case of n ao =0, p=i mm, y=1.0, 1.9 and 2.5, and 
for a*=l to 0.33. Fig. 10 confirms that the influence of y is weak in the quantity 
jRiv. \i s \ y. °> - V*(p- n/. v. n fl ) * For known n s \ the function B(a*) allows for the determi- 
nation of a relative absorption change A^i a =^ a - n ao . Fig. 10 illustrates the case ^i ao =0, but 
any other value of ji a0 is possible. The interesting point is that B(a') does not depend on 
the phase function. 
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CLAIMS 

I. A method for local and superficial (on the order of one transport mean free path) charac- 
terization of a turbid media using the following parameters: 

1) the refractive index n of the turbid media 

2) the absorption coefficient ,u a of the turbid media. 

3) the reduced scattering coefficient i± s * of the turbid media 

4) the phase function parameter*/ =( l-g 2 )/(l-gi) of the turbid media, where g { and g 2 are 
the first two moments of the phase function. 

10 and comprising the steps of: 

-measuring the spatially-resolved reflectance R(p) or any quantity which allows for the 

indirect determination of the said spatially-resolved reflectance R(p). 

-mathematically processing R(p) to compute at least one of the said parameters (n. ji^ t u s \ 

y) or the relative variations of these said parameters. 

15 -.The method of claim 1, wherein said spatially resolved reflectance is measured by an 

optical fiber probe as described by the first embodiment in section 3, and described on fig 
5 . a . 

20 
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3. The method of claim 1 and 2, wherein said spatially resolved reflectance is measured by 
an optical fiber probe composed of n fibers, as described in the particular embodiment in 
section 3, and described on fig 5.b. giving the wavelength dependance of at least one of the 
parameters (n, |a a |u s ' , y) 

4. The method of claim 1, wherein said spatially resolved reflectance is measured by a 
device put directly in contact to the turbid media, as described by the second embodiment 
in section 3 and described on fig 5.c. 

5. The method of claim 1, wherein said spatially resolved reflectance is measured by non 
contact detection unit, as described by the third embodiment in section 3 and described on 
fig 5.d. 

6. The method of claims 1 to 5, wherein the following processing is performed: 
- fit of the measured reflectance R(p) to the function: 

m,p *exp(m 3 p) (17) 

This fit gives the values for the parameters m 1? m 2 and m 3 , assuming that the expression 
r = //i ( p m: exp(/n 3 p) gives a smooth function R(p). 

The slopes ^//?(p) and ^(ln/?(p» , or any mathematical combinations of these two latter 
quantities and R(p), can be obtained directly from analytical functions using the parame- 
ters m h m 2 m 3 , or by numerical procedures from the function r = ^p^expCm.p) . 

7. The method of claims 1 to 6, wherein the absorption coefficient ji^ the reduced scatter- 
ing coefficient |i s ' and the phase function parameter y are determined on the basis of a fit 
of the measured spatially-resolved reflectance R(y, \i s \ ji at p) to a set of Monte Carlo sim- 
ulations, or functions approximating Monte Carlo simulations. 
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8 . The method of claims 1 to 6, wherein the absorption coefficient |a x the reduced scatter- 
ing coefficient |a s ' and the phase function parameter y are determined on the basis of the 
following form of the reflectance: 

R(p) = (A(p, Y, n s *) +B(^ a , n s ')) : (18) 
5 where the function A(p, y, |i s ') and B(n a , |a s ') depends also on the sources and detectors 
characteristics, and the refractive index of the sample. 

9. The method of claims 1 to 6, wherein the reduced scattering coefficient ja s * and the 
phase function parameter y are determined on the basis of the quantity J- JScpTiTTy) . 
which depends weakly on the absorption coefficient ji a for 0.3<p(i s '<5. 
10 10. The method of claims 1 to 6, wherein the absorption coefficient ji a is determined 
based on the equation: 

n a = /*[VR(p)-/(y,h s ')] (19) 
where / and k are functions that approximate Monte Carlo simulations. 

11. The method of claims 1 to 6, wherein variation of the absorption coefficient ,u a from a 

15 

known value |i ao is computed from the variation of the parameters 
jR(f>, \x s \ y, \x ao ) - jR(Q f n s \ y. v a ) • 

12. The method of claims 1 to 10, wherein said turbid medium is a biological medium. 

13. The method of claims 1 to 11, wherein the measurement an processing is repeated at 
different location of the sample (multi-site measurements), to build images of one of the 

20 

said parameters (n, ji-^ \i s \ y). 

14. An apparatus for local and superficial characterization of a turbid media comprising a 
source, a detection unit, reference means, signal processing means, at least one source 
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connected to said source unit and at least one detector connected to said detection unit, the 
distance between the source and the detector being less than one transport mean free path. 

15. The apparatues of claim 14 charaterized by the fact that the distance between the 
source and the detector varies from 0.1 to 2mm. for application to biological media and to 
turbid media similar to biological media. 

16. The apparatus of claim 14 or 15 using the method of any claim 1 to 13. 

17. The apparatus of claim 16 where the control of the homogeneity of the probed sample 
area is performed according to the procedure described in section 5.1. 

18. The apparatus of claim 16 where the calibration and normalization procedures of sec- 
tion 4. are carried out. 
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